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Twisted and rope-like assemblies of filamentous molecules are common and vital structural ele- 
ments in cells and tissue of living organisms. We study the intrinsic frustration occurring in these 
materials between the two-dimensional organization of filaments in cross section and out-of-plane 
interfilament twist in bundles. Using non-linear continuum elasticity theory of columnar materials, 
we study the favorable coupling of twist-induced stresses to the presence of edge dislocations in the 
lattice packing of bundles, which leads to a restructuring of the ground- state order of these materials 
at intermediate twist. The stability of dislocations increases as both the degree of twist and lateral 
bundle size grow. We show that in ground states of large bundles, multiple dislocations pile up into 
linear arrays, radial grain boundaries, whose number and length grows with bundle twist, giving 
rise to a rich class of "polycrystalline" packings. 

PACS numbers: 



I. INTRODUCTION 

Topological defects populate the ground states of many 
frustrated systems in condensed matter physics 0, Q • A 
key example occurs when crystalline order forms on two- 
dimensional surfaces of non-zero curvature, where the in- 
compatibility between globally straight and parallel di- 
rections generates geometrically-induced stresses that fa- 
vor defects in the crystalline order. On spheres, where 
topology requires a minimal number of twelve 5-fold 
disclinations, the problem of determining the ideal struc- 
ture of in-plane order is known alternately as the Thom- 
son or Thames problem 0,0] • This problem has impor- 
tant connections to the structure of viral capsids, which 
are closed-shell assemblies of proteins 0, ((J, and more 
recently has been the subject of experimental interest in 
the context of particle-stabilized emulsion droplets @, [1] ■ 
A clear physical picture of the coupling between sur- 
face curvature and the presence of topological defects 
has emerged based on the continuum elasticity theory 
of crystalline membranes [ii-fl2j]. In this theory, disclina- 
tions carrying a discrete topological charge act as point 
sources for in-plane stress that can be screened by a more 
homogeneous distribution of "topological charge" gener- 
ated by the appropriate Gaussian curvature of the mem- 
brane. From this viewpoint, the net topological charge of 
disclinations in the lowest energy states of curved mem- 
branes is rationally expected to increase with the inte- 
grated Gaussian curvature of a membrane, a prediction 
which has recently been tested experimentally for 2D 
crystals on surfaces of both positive and negative cur- 
vature [l3j . 

Recently, we have shown that frustration of crystalline 
order on spherically-curved surface is fundamentally con- 
nected to frustrated order in a distinct class of two- 



dimensionally ordered materials, namely, twisted fila- 
ment assemblies [H], [H| • Twisted assemblies of fibrous 
proteins are common and important structural elements 
in many biological materials, such collagen [ll| [TtJ and 
fibrin In these assemblies, helical twist of the assem- 
bly derives from the nature of interactions between chi- 
ral biofilaments, while the dense in-plane packing results 
from strong cohesive interactions between filaments 
|2~4T ] . Unlike the case of crystalline membranes where frus- 
tration arises from out-of-plane deflections, the frustra- 
tion of cross-sectional order in the twisted-filament bun- 
dle derives from a unique geometrical coupling of in- 
plane strains and filament tilts Despite the dis- 
tinct geometrical origin, twist generates stresses in the 
cross-section of filament bundles that are formally equiv- 
alent to those induced by a positive Gaussian curvature 
in a membrane, corresponding to a spherical geometry 
of effective radius, R e ff = ^ 1 /V^, where 2tt/Q is the 
pitch of helical bundle twist. Strictly speaking, due to the 
free surface at the boundary of the bundle, the twisted- 
filament packing maps more closely onto the problem 
of crystalline order of a partial, spherical cap, a prob- 
lem that has be studied theoretically for both the cases 
with (25j and without [26| topological defects. Based on 
this connection, in previous work [Til flBI ] it was argued 
that the ground-state order of filament bundles becomes 
unstable to one or more 5-fold disclinations in the cross 
sectional order, when the twist is greater than a criti- 
cal value, |0-R| C = ~ 0.47, where R is the bundle 
radius. Therefore, a range of multi-disclination ground 
states were predicted for sufficiently large and twisted 
bundles. 

In this paper, we study the continuum elasticity theory 
of twisted bundle cross sections to explore a fundamen- 
tally distinct class of topological defect configurations in 



the ground states: "neutral" 5-7 disclination pairs, or 
edge dislocations in the cross-sectional packing. Though 
the positive topological charge of a bare 5-fold disclina- 
tion best neutralizes the negative effective charge gen- 
erated by the twist, we find a broad range of conditions 
for which configurations of 5-7 dipoles, or dislocations, in 
the bundle cross section provide a lower-energy means of 
screening the geometrically-induced costs of twist. Inter- 
estingly, we find that appropriately polarized dislocations 
are universally attracted to a radial position at 
away from the bundle center, a point associated with van- 
ishing azimuthal stresses in twisted bundles. In untwisted 
bundles, as in a flat 2D crystal, edge dislocations incur an 
elastic cost per unit volume of roughly -K~o|b| 2 ln(-R/|b| ) , 
where b is the Burger's vector, Kq is the 2D Young's 
modulus [27|- In twisted bundles, we show that an elastic 
coupling between geometrically-induced stress and dislo- 
cations leads to an additional elastic energy gain per unit 
volume for optimaly-placed dislocations proportional to 
— K Q 2 R\h\. Thus, we show that twisted bundles become 
unstable to edge dislocations for reduced bundle twists 
above a threshold value, |f2i?| > (f2i?)» ~ a/Rhx(R/a), 
where a |b| is the lattice spacing of the bundle. Unlike 
the case of disclination stability [15|, find that the sta- 
bility of dislocations is governed by both the twist of the 
bundle, (£IR), as well as the size of the bundle relative 
to the microscopic size of filaments, R/a. Importantly, 
this analysis shows that the critical twist for dislocation 
stability, (flR)*, decreases, albeit slowly, to zero as the 
bundle grows macroscopic in size, as R/a —> oo. A key 
consequence of this analysis is that large bundles are un- 
stable to dislocations over a range of intermediate twist 
before 5-fold disclinations are favored for \QR\ > (Q,R) C . 

Within this intermediate range of twist (Qi£)* < 
< (flR) c , we predict a range of complex 
multi-dislocation ground states, which are quite dis- 
tinct in structure from well-separated defects in multi- 
disclination packings of twisted bundles studied in ref. 
[lij ]. In these ground states multiple 5-7 dipoles form ex- 
tended, linear chains, or grain boundaries, aligned along 
the radial directions. Grain boundaries in twisted bun- 
dles run from the free surface of bundles and termi- 
nate before reaching the central core. An example of a 
two grain boundary packing of twisted bundles is shown 
in Figure [T] The multi-dislocation ground states of 
twisted bundles are quit e similar to the "grain boundary 
scars" predicted [TTI. l28j and observed [a HH on spheri- 
cal crystals in which neutral 5-7 pairs decorate the twelve 
topologically-required 5-fold disclinations on the sphere. 
Despite the similar affinity of twist-induced stresses in 
bundles for positively-charged disclinations, at interme- 
diate levels of twist, the energetically favored defects ar- 
rays are "neutral", possessing no excess of 5-fold discli- 
nation charge. Based on numerical and scaling analysis 
of grain boundary screening of twist induced stress we 
deduce the dependence of the total number of disloca- 
tions, Nd, and number of grain boundaries, M, on twist 
and bundle size. Far above the critical twist we find that 
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FIG. 1: Cross-sectional view of the microscopic model of a 
helical filament bundle with two small-angle grain boundary 
arms, each with four dislocations. Dislocations are polarized 
such that 5-fold dislinations (red) are closer to the bundle 
center than 7- fold (blue). 

N d ~ (R/a)(nR) 2 and M ~ N d . 

The organization of this article is as follows. In Sec. 
HTl we review the continuum theory of two-dimensionally 
ordered filament arrays. In Sec. IIIII we derive the effec- 
tive energy of dislocations in the cross sections of twisted 
bundles, and determine the stability of neutral (dislo- 
cations) and charged (disclinations) topological defects 
in the ground states of twisted bundles. In Sec. HVl 
we analyze the structure and thermodynamics of multi- 
dislocation ground states and appearance of multiple 
grain boundaries in the cross sections. We also exploit a 
scaling argument based on the geometry of small-angle 
grain boundaries to establish the quantitative connection 
between the length and number of defect arrays in the 
"polycrystalline" ground states of twisted bundles. 



II. CONTINUUM THEORY OF FILAMENT 
BUNDLES 

The derivation of the equations of mechanical equi- 
librium in the non-linear continuum elasticity theory of 
twisted filament bundles have been presented in detail 
previously 15]. In this section we briefly review the key 
elements of this analysis, which allow us to construct the 
effective theory of dislocations in the next section. 

We consider a cylindrical bundle of radius, R, of fila- 
ments of unlimited length. The stress-free reference state 
is a hexagonal packing of the straight filaments in the 
cross section. We describe the elastic cost of deforma- 
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tions of the cross section order by the following energy, 



E = \ I d y{^u 2 kk + 2im,, a,,}. 



(1) 



Here, A and \x are the Lame elastic coefficients charac- 
terizing the elastic properties of the material and cor- 
respond to compressive and shear distortion of the array 
respectively, and Uij is the 2D strain of cross sectional or- 
der, defined below. In this theory, the elastic energy will 
penalize distortions of the array that change distances 
between filaments in a plane locally perpendicular to the 
filament tangent direction, t. Hence, eq. (fT]) describes 
the elastic response of hexagonal-columnar material [29j j . 
While the filaments are uniformly aligned along the z di- 
rection in the initial configuration, displacements of the 
array in general lead to filaments that are tilted into the 
xy plane of initially perfect hexagonal order 30]. This 
relationship is captured by introducing a two-component 
displacement field, u(x) = r±(x) — Xj_, describing the 
local deviation in the xy plane of a filament initially at 
x and displaced to a position r(x). For small strains, 
the in-plane displacement is related to the filament tilt 
geometrically by, t(x) = z + d z u. 

To satisfy the symmetry considerations described 
above, the non-linear strain tensor has the form 

= -(diUj + djUi + diU ■ dju - d z Uid z Ui), (2) 

which, like the displacement, only has components in the 
xy plane. The first two terms on the right hand side are 
the standard symmetric derivatives in the elastic strain 
tensor. Additionally, there are two non-linear contribu- 
tions to the strain tensor. The third term ensures the ro- 
tational invariance of the 2D solid around the z axis. The 
final term is unique to the theory of columnar materials 
and preserves the invariance of the elastic energy about 
an axis in xy plane [30j | . Since tj_ — d z u, intuitively this 
contribution to tiy shows that for a fixed separation in 
the xy plane, when neighboring filaments are tilted with 
respect to each other, the distance of closest approach be- 
tween them is reduced. The presence of this non-linear 
coupling between filament tilts and in-plane strain neces- 
sarily introduces stress in twisted filament bundles [2~i| . 

In this study, we consider helically-twisted filament 
bundles in which the cross-sectional positions of filaments 
are reorganized due to the presence of geometrically- 
induced stresses. Formally, we compose the displace- 
ment field of the helically twisted assembly from two 
deformations: an initial (^-invariant) in-plane displace- 
ment, u(x), followed by a uniform helical twist, at a 
rate Q, around the z axis. We denote the composite 
(z-dependent) displacement as uq, which has the form 

Ufi(x) = cos(fiz) [(a; + u x )x + (y + u y )y] 

- sin(fiz) [(y + u„)x - (x + u x )y] - x_i_. (3) 

In this configuration is it is straightforward to show that 



in-plane components of filament orientation have the fol- 
lowing texture, 



d z u = £lpcf>, 



(4) 



where p is the radial distance of the filament from the 
bundle center in the deformed state, and <f> is the az- 
imuthal direction, also defined with respect to the de- 
formed, or "current" position of the filament. 

The helical symmetry of these configurations allow us 
to describe the state of strain for all z, based on the in- 
plane displacement field, u(x) at z = 0. We assume the 
rate of twist and filament orientation, described by dis- 
placement uq, eq. (3), to be fixed and allow for mechan- 
ical equilibrium by relaxing the in-plane displacements 
u. This assumption has the advantage that it reduces 
the problem energy minimization to one of 2D elasticity 
theory. Minimization of the elastic energy, eq. (|T|), with 
respect to variations in u(x) leads to the Euler-Lagrange 
equations that describe the static mechanical equilibrium 
of the system, 



-djUi 



= 



(5) 



where L is the length of the bundle and the stress tensor 
has the standard form an isotropic, 2D elastic medium, 
cry = Xu^kSij + 2fiUij. Here, as in ref. [l5| we have 
neglected a term tjajkdtk/drj from eq. (|5|) because it 
contributes to the stress balance of twisted bundles at 
higher order in reduced twist, flR, which is assumed to 
be smaller than unity. As the surface of the bundle is 
free to move, we solve for states of mechanical equilib- 
rium subject to a vanish normal stress at the boundary 
of bundle 



hcrij(p = R) = 0. 



(6) 



We proceed to solve for the divergence-free stress, in 
terms of the Airy stress function x HH, related to the 
stress tensor by, 



(7) 



While this definition of cry satisfies eq. (O by construc- 
tion, it is necessary to enforce extra conditions on \ that 
ensure that the stress cor resp onds to the physical config- 
uration of u. As in ref. [10], this compatibility relation 
may be derived be equating the anti-symmetric deriva- 
tives of strain, eikejidkd^Uij , from which we derive 



s(x) + V_l x b(x) - K T . 



(8) 



where Ko = 4p(X + /i)/(A + 2p) is the 2D Youngs mod- 
ulus. The right-hand side of eq. © may be viewed as 
sources for Airy stress. The first and second of these de- 
note the sources of stress generated by topological defects 
in a bundle cross section, disclinations and dislocations 
respectively, for which the solution for u(x) is not single- 
valued. The final term, denoted as the intrinsic twist, 
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derives from the non-linear contribution to utj from fila- 
ment tilt, 



where V is the bundle volume and 



Kt = -tiktjidkdtUtj 



sn 2 



(9) 



A similar contribution derives from the non-linear cou- 
pling of elastic strain and membrane tilt in the contin- 
uum theory of elastic membranes, where it is known that 
minus the Gaussian curvature of the membrane acts a 
source for Airy stress [10(. Hence, we see that twist in 
filament bundles generates in-plane stresses that are for- 
mally equivalent to those generated by spherical mem- 
branes of curvature 3f2 2 . 

In general, two types of defects contribute to the 
right-hand side of equation © as the sources of the 
stress: disclinations and dislocations (ioj . Disclinations 
are disruptions of the orientational symmetry of the lat- 
tice and are associated with singular configurations of 
6*6 (x) = ^eijdiUj, the bond angle of the lattice. Around 
a single disclination, increases or decreases by an inte- 
ger multiple of 27r/6, 



dl ■ Vj_6» 6 = s, 



(10) 



where s = (2ir/6)n is the topological charge of the discli- 
nation. Dislocations are associated with singular config- 
urations of u(x) and defined in terms of a closed loop 
integral around which u changes by an integer multiple 
of the lattice spacing along one of the six- fold directions, 



1 • ViM t = bi, 



(11) 



where b is the Burgers vector. Multiple point defects 
in the cross section correspond to the defect densities, 
s(x) = £a««* (2) (x-xa) and b(x) = £ Q b Q( 5( 2 )(x-x Q ). 

In the presence of twist- and defect-induced stresses, 
eq. ([8]) may be solved for x and subsequently the elastic 
energy may be computed from, 



E = 



2K 



(12) 



In ref. |l5[ these equations were solved in the presence 
of an arbitrary array of disclinations in the cross section 
of filament bundles by multi-pole expansion, yielding an 
effective energy written purely in terms of charge and 
position of disclinations and bundle twist 
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- s aVi n t(x a ,Xp)sp, (13) 



V int (x a ,xp) 
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a/3 1 



16tt 2 R 2 



In 



16tt 2 V Ri)\ R? 

|Ax a/3 | 2 

(JP - pi ){R? - pl)/R? + |Ax a/3 |s 



(14) 



and Ax Q( 3 = x Q — x^. This energy has three contribu- 
tions: the first term describes the elastic cost of twist; the 
second term is the defect self energy and twist-defect in- 
teraction; and the third term describes the elastic interac- 
tion between disclinations. Importantly, both the discli- 
nation self-energy terms in (| X3[) and interaction terms 
in (|14p vanish continuously as disclinations approach the 
bundle surface, p a — > R. As noted in [l5j], this property 
derives from the screening of far-field stresses induced by 
topological defects by boundary-induced stresses as de- 
fects draw near to the free boundary. 



III. ELASTIC ENERGY OF DISLOCATIONS IN 
TWISTED BUNDLES 

A. Dislocation energies and interactions 

In this section, we take advantage of the dual descrip- 
tion of dislocations, which may be constructed from neu- 
tral 5-7 pairs of disclinations [27[ to derive the contin- 
uum theory of dislocation energies and interactions in 
twisted bundles. The theory of edge dislocations in the 
cross section of (untwisted) cylindrical crystals was orig- 
inally studied in detail by Koehler [32j. In this study, 
the resulting forms for dislocation self-energy and in- 
teraction energies were derived in terms complex area- 
integrals of stress distribution overlap, which were then 
analyzed numerically. In the present study, the exact, 
closed-form expressions for disclinations energies derived 
in ref. [l5[ allow us to derive the algebraic formula for 
the full position- and orientation dependence of disloca- 
tion energies in cylindrical crystals. 

From eq. ([5]) we may show that far field stresses gen- 
erated by a single dislocation of Burgers vector b at x, 
may be constructed by superposing a 5-fold disclination, 
s = +27r/6, at x+a/2 and 7-fold disclination, s = —2w/6, 
at x — a/2, where z x a = (2?r/6)b. Denning x+( x ) as 
the Airy stress generated by a single 5-fold disclination, 
s = +2ir/Q, which was calculated exactly in ref. [15[, the 
Airy stress corresponding to a dislocation at x, denoted 
by Xdisi(x), is given by 



Xdisi(x) = |b| lim 



X+(x + a/2)- X +(x-a/2) 



(15) 



To calculate the energy of a single edge dislocation in the 
bundle cross section, we simply superpose a 5-7 disclina- 
tion pair separated by a vector a, sum the self- and in- 
teraction energies described in eq. (1131) . and expand the 
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resulting energy to second order in a/R. This results in 
the following energy for a single dislocation, 



Edisl — E se lf + Et 



(16) 



where. 



E 



self 



VK ( 



8ir 2 R 2 



|b| 2 

8n 2 R 2 



In 



1-*-) 
R 2 ) 



a 



and 



VK 



16tt 



(17) 



(18) 



-E se z/ is the elastic energy of a single dislocation in an 
untwisted bundle, which depends largely on the radial 
position, p, and weakly on orientation of the dislocation. 
This energy is maximal for a central dislocation, p = 0, 
and reduces to the well known logarithmically divergent 
cost for a single dislocation in a bulk crystal |3^,|3^. The 
radial dependence of E se if shown in Fig. [2^, becomes 
singular as the dislocation approaches the bundles surface 
as the boundary-induced force on a dislocation diverges 
as ~ (R — p) _1 [HJ- Hence, in the limit that R— p <C a, 
the 5-7 disclination superposition is non-analytic asa-> 
0, and hence the small-a expansion of eq. ([T71) becomes 
inaccurate. 

Etwist describes the elastic coupling between the twist 
induced stresses and the dislocation stresses shown in 
Fig. [2b. Notably, this coupling is negative and minimal 
for b = b<f) and p — R/y/3, demonstrating that twist 
favors dislocations of a certain polarization and located 
at a specific radial position in twisted bundles. The fa- 
vorable orientations of dislocations correspond to discli- 
nation dipoles oriented along the radial direction, with 
the 5-fold end oriented towards the bundle center. Alter- 
natively, we may view such a dislocation as a partially 
removed row of filaments extending from the free bundle 
surface to the dislocation. 

To understand the origin of an optimal location of dis- 
location in twisted bundles, we consider an alternative 
derivation of E twist based on the Peach-Koehler force [35[ 
generated by twist-induced stress. The force per unit 
length on a dislocation line along z subject to imposed 
stress (Tij is given by, 



PK 



(19) 



The stresses generated by twist are described by the so- 
lution to Kq V\xtwist = —3Q 2 , which can be readily 
solved to show the following azimuthal stress distribu- 
tion, 



rl7 St ( P ) = ^-(R 2 -5p 2 ). 



(20) 



This stress distribution divides the bundle into two re- 
gions: tensile hoop stresses, a l V? st > 0, at the bundle core 
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FIG. 2: The dislocation self-energy is plotted in (a) vs. radial 
position, p. In (b), we show the radial dependence of the elas- 
tic coupling between twist and dislocation stresses. In (c), a 
schematic showing the Peach-Koehler force on disclocations 
in the presence of twisted induced stress, where _L indicates 
the position of an edge dislocation. The dashed line indicates 
a contour of vanishing hoop stress, to which dislocations in 
highly-twisted bundles are driven. For each figure, the dislo- 
cation orientation is b = bd>. 



for p < Rj -\/3; and due to large azimuthal tilt of filaments 
at the periphery, compressive hoop stresses, a 1 ™™ 1 < 0, 

for p > R/y/3. Since twist induces a radially symmetric 
stress, cr^ lst = 0, the force of a dislocation whose Burg- 
ers vector is oriented along <fi (with a 5-7 dipole along 
f) is in the r direction. As shown in Fig for such a 
defect in the compressive zone, for p > i?/v3, the Peach- 
Koehler force drives the dislocation inwards, while in the 
tensile zone, for p < R/y/3, this force drive the dislo- 
cation outwards. Hence, the force vanishes where o 1 ^™ 1 

vanishes at p — R/y/3, the stable position. Combining 
eqs. (TT9"|) and (|2Tj)) . energetic coupling between twist and 
dislocations, eq. (IT51) . may be readily calculated from 
the mechanical work of driving a defect into the bundle, 



We end this section with an analysis of dislocation- 
dislocation interactions in cylindrical bundles. As with 
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E se if and E tW i S t, we derive these from the interactions 
between two neutral disclination pairs. To compute inter- 
dislocation energies, we sum disclination interactions 
over two disclinations, sf — ±27r/6 at xf = x± ai/2, 



and the second pair of disclinations, sjf = ±27r/6 at 
x^ = x ± a 2 /2. Again, we retain terms to lowest order 
in a from the expansion of multiple disclination interac- 
tions, yielding the interaction energy 



Ej. 



1 



K V A-k 2 R 2 



, w, 9 * 9 ^\ (ri x bi)(ro x bo) a ^ 
(bi • b 2 ) (In cos 2 £ + sin 2 £) + ^ ^ sin 4 £ 



(bi x Axi 2 )(b 2 x Ax 2 i) 



R 2 

(bi x Ax 12 )(b 2 x r 2 )(l- p\/R 2 



|Ax 12 p 

(b 2 x Ax 21 )(b 1 xr 1 )(l-p 2 /i? 2 ) 



(1 



(il 2 - pf)(E 2 - p|) + I Axi 2 | 



sin 2 £ 



(21) 



Here, r$ measures the position of ith dislocation with 
respect to the bundle center, and £ is defined by 



cos 2 £ 



|Ax 12 | 



(i? 2 -p 2 )(i? 2 -p 2 ) + |Ax 12 | 2 ' 



(22) 



Due to the presence of the free boundary, this pair poten- 
tial encodes a significantly more complex dependence on 
defect orientation and position than the well-known elas- 
tic interactions of dislocations in 2D crystals [27], |33|, 134J . 
However, we notice the well-known form of logarithmic 
dislocation interactions in bulk crystals is easily obtained 
by taking the limit that pt/R — > of eq. (PHI) , for 
which cos£ — !• |Axi 2 |/i? and sin£ — > 1. Additionally, 
we note that when in the limit where either dislocation 
approaches the boundary, pi/R^ 1, the dislocations in- 
teractions vanish, which can easily be verified for the case 
cos£ — > 1 and sin£ — > 0. 




FIG. 3: The phase diagram indicating stability of appropri- 
ately oriented (b = b(j>) dislocations and 5-fold disclinations 
in cross sectional order of twisted filament bundles in terms 
of reduced twist and reduced size of bundles. 



B. Defect phase diagram of twisted bundles 

Here, we analyze the stability of disclination and dislo- 
cations in the cross section of twisted bundles. As shown 
previously [3, HU, and eq. (fT3|) . twist-induced stresses 
couple favorably to the presence of positively charged 
(s = +2"7r/6), 5-fold disclinations, and above a critical 
threshold of reduced twist, |(fi-R) c = y/2/9 — 0.47, this 
energetic coupling is sufficient to compensate for the pos- 
itive self-energy cost of a single disclination at any posi- 
tion. Thus, bundles are unstable to one or more 5-fold 
disclinations for \SIR\ > (SIR) C . 

We consider the stability of a dislocation by consider- 
ing the energy of a single dislocation, which is polarized 
by twist-induced stress such that b = b<j>. Minimizing 
the sum of eqs. (fTTl) and (fT8)l over radial position, p, we 
find the value of twist, (SIR)* at which the net cost of a 
single dislocation vanishes, Ediaiip*) = 0, where p* is the 
stable position of the dislocation [43| For larger bundle 
twists |f2i2| > (fiii)*, one or more dislocations is stable 
in the low-energy packing of twisted bundles. 

In Fig. [3] we show the value of both the threshold for 
disclinations and dislocations, (CIR) C and (Oi?)*, respec- 



tively, as functions of R/a, the size of the bundle in units 
of the lattice spacing, a ~ |b|. The threshold for 5-fold 
disclinations is independent of bundle size; however, we 
find that (SIR)* < (QR) C for all R/a > 2. Thus, for fixed 
bundle size R/a, for increasing values of (SIR), twisted 
bundles become unstable to neutral defects, dislocations, 
before becoming unstable to the 5-fold disclinations in 
their ground-state packing. 

We can roughly estimate the size dependence of (SIR)* 
in the regime of large bundles. In this limit, the position 
of the dislocation is determined by the twist energy alone, 
which is minimal for p* — R/y/3. Solving Edisi(R/ VS) = 
critical twist (SIR)*, we find 



irR 



In 



R 



0.072 



(23) 



This formula highlights the balance between the logarith- 
mic self-energy of a single dislocation, ~ ifo|b| 2 ln( R/a), 
and the compensating dislocation-twist coupling, ~ 
—KQSl 2 R\h\. Hence, we find that the threshold twist 
necessary for stabilizing dislocations in the cross section 
becomes arbitrarily small as bundles become macroscopic 
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in radius, in the R/a —> oo limit. This analysis suggests 
that dislocations proliferate in large bundles, and there- 
fore, understanding the ground state packing requires the 
study of multi-dislocation structure and energetics. 



IV. MULTI-DISLOCATION GROUND STATES 
A. Numerical study 

In this section we explore the structure of multi- 
dislocation cross sections in the region of intermediate 
twist, (f2i?)» < \CIR\ < (£IR) C . We base our analy- 
sis on a certain class of mechanically stable and high- 
symmetry dislocation geometries where parallel, b = bcj) 
dislocations concentrate along M identical radial lines, 
or "arms", spaced evenly at angular intervals of 2tt/M, 
around the bundle. A similar class multi-dislocation 
geometries have been studied in the context of grain- 
boundary screening of isolated disclinations in 2D crys- 
tals (36|. In these geometries, each radial line of dislo- 
cations is line of mirror symmetry in the defect packing 
so that 0>0 = along these lines and, by eq. (Q33), the 
(j> component of force (the glide direction) vanishes for 
each dislocation. The remaining force balance along the 
radial direction results from repulsive inter-dislocation 
forces that favor expansion of the array and the Peach- 
Koehler force on dislocations arising from twist-induced 
stresses that favors a restoring compression of the array. 
As shown in Fig. [TJ extended strings of alternating 5- 
and 7-fold defects constitute tilt grain boundaries across 
which the orientation of two domains of crystalline order 
rotates by a discrete angle [34| . 

To determine the radial position of dislocations in these 
minimal-energy configurations, we consider the total en- 
ergy of configurations possessing total disclinations, 
composed of M equivalent arms of n = Nd/M disloca- 
tions per arm. For a fixed dislocation geometry, reduced 
twist and bundle size, the sum of the single defect energy, 
cq. (|16p. and interaction energy between defect pairs, eq. 
(1211) , is numerically minimized with respect to the radial 
position of the n dislocation "rings" in the array. In this 
analysis, the minimum spacing between successive dislo- 
cations along the array is set to be, a, the lattice spacing. 

Fig. H] shows the results for the number and arrange- 
ment of dislocations in a bundle of size R = 100a for 
a range of twist below the threshold for stable 5-fold 
defects. As the twist increases beyond (QR)„ ~ 0.16, 
the number of dislocations favored in the cross section 
increases quickly. For each value of Nd, the geometry 
of the dislocation packing is labeled by the integer pair, 
(M, n), denoting the number of grain boundary arms and 
the number of dislocations per arm, respectively. Along 
with the total dislocation number, the number of radial 
grain boundaries also grows with (CIR), leading 2-, 3-, 4- 
and 6-fold grain boundary geometries depicted in Fig. 2) 

In Fig. [S]we show results for Nd vs. (flR) for a much 
larger bundle, R = 700a. While this bundle shows a simi- 
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FIG. 4: (Top) The total dislocation number for ground-state 
configurations of twisted bundles with multiple dislocations 
for a bundle of size, R — 100a. Integer pairs, (M,n), refer 
to the number grain of boundary arms and the number of 
dislocations per arm, respectively. (Bottom) A, B, C, and D 
show 2-, 3-, 4- and 6-fold grain boundary geometries, where _L 
labels the position of a single dislocation. 



lar trend with increasing twist, we note that the threshold 
for stable dislocations is markedly reduced, (fii?)* ~ 0.07 
and a distinct sequence of dislocation geometries is pre- 
dicted as Nd increases rapidly with twist. Notably, we 
find for all multi-dislocation geometries over a range from 
R/a = 20 to 700, that grain-boundary arms penetrate 
only a fraction of the distance from the bundle surface 
to the bundle center, terminating in the bulk at a finite 
radius, a feature uncommon in bulk crystalline materials. 

To investigate the evolution of grain-boundary struc- 
ture in bundles with increasing twist, in Fig. [5] we plot 
the number of grain boundary arms, M , vs. Nd for all 
values of R/a studied. Over the range of dislocations 
explored here (up to Nd — 50) we find little systematic 
dependence of the growth in the number of grain bound- 
aries on bundle size. Despite more than an order of mag- 
nitude variation of bundle size, the trend of increasing 
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FIG. 5: The total dislocation number for ground-state con- 
figurations of twisted bundles with multiple dislocations for a 
bundle of size, R = 700a. Data points are labeled as in Fig. 
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FIG. 6: Plot of number of grain boundarm arms, M, vs. 
total dislocation number, Nd for bundle sizes in the range 
R/a = 20 - 700. 



number of grain boundaries is consistent with a roughly 
linear relationship, M ~ Nd for all R/a. This suggests 
that the optimal dislocation geometry is nominally de- 
termined by Nd alone, which in turn is regulated by flR 
in a manifestly size-dependent manner, as evidenced by 
the results for R = 100a and R = 700a bundles, of Figs. 
Hand El respectively. 



B. Scaling analysis 

Above, we found that the gross structure of the multi- 
dislocation ground state is predominantly sensitive to 
the total dislocation number. Here, we consider a sim- 
ple scaling argument to understand the dependence of 
Nd on bundle twist and size. This argument is similar 
to the geometric analysis of "grain-boundary scars" on 
spherical crystals [HI], with the notable exception that 
in the present case, neutral grain boundaries form in the 



absence of excess point disclinations. According to the 
compatibility relation, eq. (jSJ), we can formally consider 
the source of twist-induced stresses, Kt = 3£1 2 , to be a 
uniform areal density of negatively charged disclinations. 
Integrating this charge density over the cross section of 
the bundle we define an effective disclination charge, 



s eff = -3ir{QR) 2 



(24) 



As argued in ref . [llj , the strain generated by this "topo- 
logical defect" can be compensated by the presence of M 
radial grain boundaries, each of which accommodates a 
rotation of 9 ~ a/d, where d is the mean dislocation 
spacing along the boundary (3?]]. Equating the effective 
topological charge to the total grain boundary rotation, 
we find the mean-spacing between dislocations, 



or 1 (sir) 2 /m. 



(25) 



Integrating the linear density of dislocations along the 
length of grain boundaries (~ R) we find the mean num- 
ber of dislocations per arm, 



n « (R/a)(QR) 2 /M. 



(26) 



Multiplying n by the number of grain boundaries in 
the cross section, we argue that for large twist, Nd ~ 
(R/a)(SlR) 2 . To capture both limiting cases of large 
twist and the critical twist at which Nd vanishes, we 
construct the following scaling form for total disclination 
number, 



N d ~(R/a)[(QR) 2 -(SIR) 2 



(27) 



Hence, not only do larger bundles become unstable to 
dislocations at smaller values of bundle twist, the growth 
of the optimal number of dislocations with "excess" twist 
in large bundles is also more rapid than in smaller bun- 
dles. 

In Fig. [7] we compare the total dislocation number of 
the numerically-determined ground states to the scaling 
prediction, eq. ([27]) . We find that the scaling prediction 
agrees well with numerical results over the range of Nd 
and the large range of bundle sizes explored here, 20 < 
R/a < 700. 



V. DISCUSSION 

In summary, we have shown that geometric frustration 
arising from helical twist in two-dimensionally ordered 
filament bundles restructures the ground-state packing 
at intermediate twist by favoring the presence of appro- 
priately oriented, edge dislocations. Based on the con- 
tinuum theory of disclinations and dislocations in fila- 
ment bundles, we show that dislocations become favor- 
able in the cross section at twist smaller than the crit- 
ical twist needed to stabilize 5-fold disclinations. Un- 
like the case of stable 5-fold disclinations studied previ- 
ously [3, Ell) here we find that the threshold for sta- 
ble dislocations in bundles is highly dependent on the 



9 



• 


. 




• It/ a — ( UU 


•sss: 
— ~ 


• R/a = 500 ' 


-ssr 


• R/a — 4UO 




• R/a = 350 




• R/a = 250 ; 




• R/a = 150 : 




• R/a = 100 '. 




■ R/a = 50 ■ 




• R/a = 20 ■ 





5 10 15 20 25 



(R/a)[(nRf - (SlR)l] 

FIG. 7: Plot of dislocation number vs. scaling prediction, eq. 
l|2?). for bundle sizes in the range R/a = 20 - 700. 

size of the bundle compared to the microscopic inter- 
filament spacing. Above the threshold single dislocations 
in twisted bundles, we predict a rich spectrum of low- 
symmetry ground-state order. Twist-induced stresses in 
filament bundles lead to a natural tendency to "polyg- 
onalize" the cross-sectional packing, giving rise to low- 
energy structures where multiple crystalline domains are 
separated by radially-extending grain boundaries that 
terminate in the bundle bulk, not unlike the finite-length 
grain-boundary scars of spherical crystals. 

This study is significant in the context of frus- 
trated order because it demonstrates that "neutral" 
configurations of topological defects may effectively 
screen "charged" , geometrically-induced stresses, like the 
stresses generated by filament twist in bundles. Previous 
studies of defects on curved, crystalline membranes, have 
predicted extended chains dislocations, or scars, only in 
the presence of excess disclinations that are themselves 
either forced in by topology [ill HH ° r j i n the case of 
a membranes with a free boundary, as the result of en- 
ergetic coupling to curvature-induced stresses [25], l38j . 
In these cases, the dislocation arrays function to screen 
the disclination stress more efficiently than the stresses 
induced by Gaussian curvature. Here, we show that neu- 
tral arrays of 5-7 disclination pairs flood the ground-state 
packing of crystalline bundles well before twist favors 
the incorporation of excess 5-fold defects. That is, dis- 
locations arrays are also driven into the packing of frus- 
trated materials by the tendency to screen geometrically- 
induced stresses alone. Due to the formal relationship 
between the non-linear elasticity of twisted bundles and 
curved, crystalline membranes, we expect that the novel 
grain-boundary geometries predicted here may also occur 
as ground states of the latter system. 

In the context of filamentous materials, the present 
study is significant for two reasons. First, it identifies 
flR and R/a as the two geometric parameters that gov- 
ern the ground state packing of helically-twisted bundles. 



Importantly, we show that the critical degree of twist at 
which the bundle cross section becomes unstable to topo- 
logical defects is crucially sensitive to bundle size, R/a. 
We may relate the reduced twist of the bundle to the tilt 
angle, 9 max of the outermost filament with respect the 
pitch axis of the helical bundle by, 9 max = tan -1 (ili?). 
For small bundles, less than a few radial filament layers, 
the critical twist for stability of any defect type in the 
cross section corresponds to a degree of tilt greater than 
25°. Because the critical twist decreases with bundle size 
as (Sli?,)* ~ a/Rln(R/a), for bundles that are macro- 
scopically large compared to filament size, say R = 100a, 
the critical degree of filament tilt is markedly reduced to 
nearly 9°. For comparison, we note the helical twist of 
certain collagen fibrils is in the range 15 — 17° of helical 

tilt [TSUI. 

The case of collagen points to the second important re- 
sult regarding the structure of optimally-packed twisted 
fibers. In many tissue types, collagen fibrils form with 
lateral dimensions hundreds of times larger than the 
roughly 1 nm scale of constituent filaments, making mul- 
tiple dislocations energetically favorable even in fibrils of 
relatively modest twist. Notably the precise nature of 
the cross-sectional ordering of collagen molecules in fib- 
rils is a long-standing and open question, in part, due to 
small-angle scattering data that suggest cross sections are 
composed of unknown superpositions of crystalline inter- 
molecular order and disordered inter-molecular packing 
of some type. Numerous models have been proposed 
to infer the real-space packing [40l ]. many of which mix 
aspects of crystalline and non-crystalline order in novel 
ways [4lj . To date, the model most consistent with ob- 
served features of x-ray scattering data was proposed by 
Hulmes, Wess, Prockop and Fratzl [42j . In this model, 
multiple crystalline domains in a cylindrical bundles are 
separated by grain boundaries extending radially from 
a central, low-density region to the surface of the fibril. 
Remarkably, this model is very similar in gross struc- 
ture to the ground states of twisted bundles predicted for 
large bundles of intermediate twist (for example in Fig. 
U). Though this model of collagen fibril packing did not 
take into account the effects of twist explicitly, we believe 
many of the key features of the "disordered" packing of 
these materials may be understood as crucial elements 
of energy-minimizing packings of twisted bundles. Fu- 
ture work will explore the form factor of ideal packings 
of twisted bundles and critically test the intriguing and 
putative connection between defects in the ground states 
of twisted bundles and the disorder in the collagen fibril. 
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